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We have studied the electronic structure of the spin-gapped 
system CSV2O5 by means of an ab initio calculation. Our 
analysis and a re-examination of the susceptibility data indi- 
cate that the behavior of this system is much closer to that of 
an alternating spin-i antiferromagnetic chain with significant 
inter-dimer coupling and weaker inter-chain couplings than 
that of isolated dimers as was initially proposed. Comparison 
to the vanadate family members a'-NaV205, 7-LiV205 and 
isostructural compounds like (VO)2P207 is discussed. 

PACS numbers: 75.30.Gw, 75.10.Jm, 78.30.-j 



I. INTRODUCTION 



Low-dimensional spin-gapped quantum systems are 
of current interest since they show interesting ground 
states and a variety of unconventional low-lying exci- 
tations. Examples of such systems indude spin-1 Hal- 
dane chainsEl, spin-i even- leg ladderso and spin-i al- 
ternating chainscl. The discovery of appropriate conu 
pounds likej^2BaNi05cl (Haldane chain), Sri4Cu2404iEl 
or SrCu203El (ladder systems) has brought new insight 
into the study of these systems. There is a long list 
of various spin-i alternating chain systems which are 
being intensively studied in connection to their mag- 
netic excitations, to mentiou. a few of them KCUCI3Q, 
TlCuClala, Cu(N03)2.|D20y. If frustration is also 
present new features appear in the magnetization spec- 
trum of|-the low-dimensional ™«intum systems— jas in 
CuGeOalia, (VO)2P207 (VOPO)^, SrGaj2(B03)2li3 or in 
the recently synthesized Cu2Te205Br2E^. 

An important family of low-dimensional compounds 
arcpthe layered vanadates AV2O5 (A^Ca, Mg, Na, Li, 
Cs)EJ. While CaV205 and MgV205 contain only mag- 
netic V''"'" ions and behaveJike spin- ijadders with spin- 
gaps of the order of GOOKB and 2OK0 respectively, a'- 
NaV2 05, 7-LiV205 and CSV2O5 are mixed- valence sys- 
tems (V^'^"*" on average) with important charge and spin 
fluctuations. This last property makes this family spe- 
cially interesting since it allows to study the influence of 
charge-ordering and the corresponding crystallographic 
distortions on the magnetic interaatioiw. The quarter- 
flUed ladder compound a'-NaV205llII't§ is a highly dis- 
cussed material and has become a model substance for 
the study of spin, charge and orbital coupling. It be- 
haves like a spin-i Heisenberg system at high tempera- 



tures and at 34K undergoes a charge ordering (2V^ ''^+ 
— > -fj_L^^+) transitionO accompanied jsath a lattice 
distortionEJ and the opening of a spin gapcJ. 7-LiV205 
is a charge ordered, system with no spin gap opening at 
low temperaturef3. Recent investigations rbased on an 
ah initio analysis of the electronic structuraH^ show that 
charge fluctuations in this systenj-are still signiflcant and 
the experimental observationg23'Ej can be well described 
by a spin-i asymmetric quarter-filled ladder model. Such 
proposal has also found support by ja_numerical analysis 
of optical conductivity experimentsEj. The third mem- 
ber of this vanadate family, CSV2O5, has been given 
much less attention— Existing susceptibility measure- 
ments of Isobe et alE3 could be rather well explained by 
considering this material as a system of isolated dimers 
with a spin-gap of « 12meV. Measurements other than 
that are unkiiown to us except for older susceptibility 
measurement£3 which somewhat disagree with Isobe et 
al.'s data and different values of the exchange integral J 
and the gyromagnetic ratio g are predicted. 

There are various reasons why we believe that this 
system deserves attention: (i) the question of how the 
electronic and magnetic properties of this material com- 
pare with the properties of the other two members of the 
alkali- vanadate family, i.e. a'-NaV205 and7-LiV205 (ii) 
this system crystallizes in a similar monoclinic structure 
as the most discussed VOPO and many other materials 
e.g. KCuCla, TICUCI3 showing an alternating chain be- 
havior of the magnetic interactions. It is then worthwhile 
to find out whether CSV2O5 has also a similar behavior 
of the magnetic interactions as its isostructural members, 
(iii) as will be shown below, a re-examination of the sus- 
ceptibility measurements reveals that a dimer model is 
not the only model that can provide reasonable fitting of 
the susceptibility data. 

In view of the above, we have performed a micro- 
scopic study of this system by considering Density- 
Functional-Theory-based (DFT) calculations. The mi- 
croscopic study of the electronic behavior of this com- 
pound reveals that the inter-dimer coupling in this ma- 
terial in addition to the intra-dimer interaction is also sig- 
nificant and shouldn't be neglected. Further, an analysis 
of the susceptibility data shows that getting estimates 
of the coupling constants and the gyromagnetic factor 
based merely on a fitting procedure is very subtle and 
with such an analysis one can only test the consistency 
with the assumed model, but cannot prove the unique- 
ness of that model. Results on that will be presented in 
the next sections. We shall also discuss the similarities 
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and differences to the other members of this vanadate 
family, Q;'-NaV205 and 7-LiV205 as well as to isostruc- 
tural compounds like VOPO. 



II. CRYSTAL STRUCTURE 

The crystal structure of CSV2O5 is somewhat different 
from that of a'-NaV205 and 7-LiV205. a'-NaV205 and 
7-LiV205 crystallize in the orthorhombic space group 
Pmmn (£'2^ ^^'^ Pnma (£'2^ respectively. They con- 
sist of a double chain-structure of edge-sharing distorted 
square VO5 pyramids running along the orthorhombic y 
axis, which are linked together via common corners to 
form layers. These layers are stacked upon each other 
along z. Na/Li atoms are located between these layers. 
CSV2O5 crystallizesEII, on the other hand, in the mono- 
clinic space group P2i/c (C|^). CSV2O5 has also a lay- 
ered structure with Cs ions in between the layers. These 
layers, which are somewhat tilted ( i.e. not strictly in 
yz plane), and are stacked upon each other along x, are 
made of two types of vanadium atoms VI and V2. The 
VI atoms are in a distorted square pyramidal environ- 
ment of oxygens and form pairs by sharing the edge of the 
pyramids. These pairs are separated from one another by 
the V2 vanadiums in a tetrahedral co-ordination of oxy- 
gens as shown in Fig. ^ Within a pair of square pyra- 
mids, the two apical O atoms are pointing in opposite di- 
rections and the remaining O atoms of the pyramids are 
shared with tetrahedra. The Vl-Vl distance across the 
shared edge of the square-pyramids-pair is 3.073 A and 
the V1-V2 distance across the shared corner of a pyra- 
mid and a tetrahedron is 3.352 A. Since the average V2- 
O distance within the tetrahedron d=1.718 A is shorter 
than the average Vl-0 distance within the square pyra- 
mid d=:|l_882 A, one can conclude from a bond-valence 
analysisEZl that the VI sites are in the oxidation V^^ and 
the V2 sites are in the oxidation V^"*". As the intra-pair 
V^^-V^^ distance is 3.073 A while the inter-pair distance 
along the z-axis is 5.501 A, this compound has been as- 
sumed to be a system of isolated dimers. 



III. BAND STRUCTURE 

In order to analyze the electronic behavior of CSV2O5, 
we carried out DFT calculations within the Lo- 
cal Density Approximation (LDA) by employing the 
fuU-poteptial linearized augmented plane wave code 
WIEN9728 (LAPW) and by hncar-mufiin-tin orbital 
(LMTO)Ei based on the Stuttgart TBLMTO-47 code. 
Both calculations agree well within the allowed error- 
bars of the various different approximations involved in 
two different methods. In both approaches we treated 
the exchange-correlaticp^ part by using the generalized 
gradient approximation^!]. 



In Fig. I the energy bands for CSV2O5 along the path 
TBDZTYE are shown. Since the unit ceU of CSV2O5 
contains four formula units, i.e. four V^^" and four V^+, 
there are four bands around the Fermi level which corre- 
spond to V^^-Sd orbitals. These bands are separated by 
an energy gap of about 2 eV from the lower 0-p valence 
bands and an energy gap of about 0.5 eV from the up- 
per conduction bands. The conduction bands up to 4eV 
are of V-3d nature. There are a few points to be noted 
here: (i) The energy range on the CSV2O5 banck-js of 
the same pwicr of magnitude as it is in Q;'-NaV205liZl and 
7-LiV205Ea. (ii) The contribution of V^+-3d is mostly 
in the upper conduction bands, starting from about 0.75 
eV above the Fermi level and have a very small contri- 
bution at the Fermi level indicating that, contrary to 
a'-NaV205 and 7-LiV205, here the system is close to 
complete charge ordering, (iii) There is no dispersion 
along the x axis while along the z axis (chain direction) 
the dispersion is non-negligible and of about 300meV. 
Along the y axis the LDA-bands show a small disper- 
sion which corresponds to inter-chain interactions. In 
order to investigate the dimension-dependent behavior 
of CSV2O5, we show in Fig. ^(a) the partial density of 
states (DOS) of the Y'^+ 3d-orbitals for CSV2O5. For 
comparison we have also plotted in Fig. |^(b) the partial 
density of states of Vl-3d orbitals for 7-LiV205E3 and 
in Fig. |(c) that of V-3d orbitals for a'-NaV205. Note 
that, while the 7-LiV205 and Q;'-NaV205 DOS show the 
characteristic quasi-lD van Hove singularities near the 
band edges around the Fermi level, this is not the case 
for the Cs-compound. The DOS for CSV2O5 shows more 
a 2D-behavior than ID. (iv) In the global frame of refer- 
ence, the character of the bands around the Fermi level 
turn out to be that of d^^ with small admixture from 
da;2„j,2 character. Rotation to local co-ordinate system 
with the z-axis aligned along the (V-apical O) direction 
for individual pyramids transform this to dxy character, 
in agreement with the crystal-field analysis of a V*^ ion 
in the pyramidal environmendU. 

The band-structure analysis leads us to the conclusion 
that inter-dimer coupling is non-negligible. If the system 
would have been a pure dimer system we would have 
expected no dispersion in any direction. The existence 
of a non-negligible dispersion tells us that a model of 
unlinked dimers is too simple and there must be some 
inter-dimer coupling to explain the band picture. From 
the DOS analysis we learn that inter-dimer interactions 
in both z and y direction are important for the electronic 
properties. The possible role of the V^+04 tetrahedra in 
this context will be investigated in the next section. 



IV. MICROSCOPIC MODEL (EFFECTIVE 
MODEL) 

In order to understand the nature of the ground state 
in CSV2O5 we need to determine the appropriate mi- 
croscopic model which explains the low-energy physics 
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of this compound. We have therefore apphed the so- 
called downfolding method offered by the new and gen- 
eralized version of the LMTOEil method, together with 
a tight-binding (TB) analysis on the band-structure re- 
sults. The downfolding technique consists of integrating 
out the high energy degrees of freedom so as to describe 
the details of LDA energy bands close to the Fermi energy 
in terms of few-orbital effective Hamiltonians. From this 
analysis, we can extract the effective hopping matrix ele- 
ments tij between vanadium ions in CSV2O5 (by Fourier 
Transform of the downfoldcd Hamiltonian H{k) — ^ H{R) 
) which reproduce the behavior of the LDA bands close 
to the Fermi level. Since, as seen in the previous sec- 
tion, the contribution of the V^^-Sd orbitals around the 
Fermi level is very small, and the biggest contribution 
comes from the Y'^^-'idxy orbital, a one-orbital effective 
tight-binding Hamiltonian should give a good description 
of the low energy physics. The tight-binding Hamiltonian 
therefore, can be written as: 

HTB^-Y.^-M'^^ + ^tc,) (1) 
<i,j> 

where i and j denote a pair of V'^ ions labeled each one 
by (n, s) where n denotes the unit cell and s — 1,2,3,4 
V*"*" ions in the unit cell. In Fig. ^all the considered hop- 
pings ti, t2, and t5 have been drawn. The subscript 
1, 2, 3 and 5 indicate respectively that the hopping in- 
tegral is between first nearest neighbors (n.n.), second 
n.n., third n.n and fifth n.n. The hopping term t^, which 
connects two V^"*" ions belonging to two different layers 
along X, has been neglected due to its vanishingly small 
value. Note here that no hopping terms have been con- 
sidered along the x direction, since as we have learnt from 
the band picture, there is no dispersion along this direc- 
tion and that means that the interlayer coupling along 
X must be negligible. In Fig. ^ we show a comparison 
of the downfolding-TB bands to the LDA bands. We 
observe that not only the hopping corresponding to the 
intra-dimer interaction ti is important but also the inter- 
dimer matrix element ^3 along z is considerable. A model 
with only ti and t^, i.e. alternating chain, shows a dis- 
persion relation of the form: 



E = ±^tl + tl + 2tit3Cosk^ (2) 

which is symmetric with respect to E=0 (i.e. to Ep) and 
doesn't account for the band-splittings in the paths T — D 
and r — Y. These features are described by the inter- 
chain hopping parameters along the y direction, ^2 and 
^5, though they are much smaller than ti and ^3. Note 
that the TB-rcsults show an artificial band-crossing along 
the path B — D and Z — T. We have analyzed this result 
and seen that the inclusion of longer-ranged interactions 
lifts this artificial crossing. 

Concerning the nature of the hopping integrals, we ob- 
serve that the values of the dominant hopping parameter 
ti — Q.12eV (which agrees with the overall band- width 



of the V^^-daiy derived low-energy bands in CSV2O5) is 
smaller titan the dominant hopping parameter for a'- 
NaV205tZl ta = 0.37ey and 7-LiV205Ea ta = 0.35ey. 
This difference can be explained by the nature of the 
path. While ti in CSV2O5 corresponds io aV — O — V 
path in between two ed(7e-sharing square pyramids, ta in 
a'-NaV205 and 7-LiV205 corresponds to a V — O — V 
path between two corner-sharing square pyramids. This 
feature has important implications in the strength of the 
exchange interaction as has been pointed out in Ref.Lj. 
Also note that the same kind of ed(7e-sharinffi-4)yramids 
are also present in a'-NaV205 and 7-LiV205E3 with ef- 
fective hopping values smaller than ta and similar to ii 
in CSV2O5. 

About the nature of the rest of the hopping integrals in 
CSV2O5, we observe that t2, ^3 and is correspond to 

paths through the bridging V^+04 tetrahedra with 
different lengths and angles. Beltran et a/.E3 studied, 
based on geometrical considerations, the various super- 
exchange paths for the case of vanadyl phosphates. In 
particular they considered the VOPO compound. VOPO 
and CSV2O5 have a similar monoclinic structure if one 
identifies the P^+04 tetrahedra and the V^+05 square 
pyramids in VOPO with the V^+04 tetrahedra and the 
^^+05 square pyramids in CSV2O5 respectively. As ar- 
gued by Beltran et ai, based on the formation of coher- 
ent molecular electron orbitals, the coupling between two 
V^+ through bridging oxygens for separated pyramids is 
weaker than those via the pathways involving V^"*" ions. 
Furthermore, double bridging modes, involving two V^"*" 
ions, as in t2 and t^ (see Fig. ^ can be stronger than sin- 
gle bridging modes involving one V^+ ion as in ts. These 
interactions can be as strong as paths of the type ti, i.e. 
V^''"-0-V*"^ where the square pyramids are linked to each 
other by sharing an edge. Geometry considerations then 
lead to a relative order of importance among the vari- 
ous bridging modes, which for CSV2O5 is observed to be 
^3 > ^5 > ^2- We note here that while in CSV2O5 we have 
seen that the dimers are formed by the edge-sharing pyra- 
mids (ti) with the predominant inter-dimer interaction 
formed by V*''"-V^"'' interaction bridged by V^+04 {t^), 
in VOPO these roles are interchanged in the sense that 
dimers are forped,by V^+Os square pyramids bridged by 
P^+04 groupaHa {t^) and the inter-dimer interaction is 
formed by neighboring V*"'"05 pyramids (ti). This is due 
to the different type of distortion in the square pyramids 
as well as different angle and distance relations between 
these two compounds. 

Using the relation Ji = At\/U ^ U being the on-site 
Coulomb electron-electron interaction and considering 
U~ 2.8 eVll3 the DFT calculation gives an estimate of 
the intra-dimer exchange coupling, Ji 225K . 

To summarize, the conclusions that we draw from the 
above analysis is that ab initio calculations support an al- 
ternating chain behavior for CSV2O5 along z with weaker 
inter-chain interactions along y. 
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V. SUSCEPTIBILITY 



VI. CONCLUSIONS 



In order to test the results obtained from the ab initio 
calculation-!^ only have available two sets of suscepti- 
bility dataE3'Ea as a function of temperature which don't 
completely agree quantitatively with, each other. We con- 
sider here the data by Isobe et al.tn These authors pro- 
posed that their experimental data are well fitted by the 
susceptibility of a spin-i dimer system with H — JS1S2 
given by: 



'X.raw 

1 



XCsV-iO;, + Ximp 
XO + Ximp 



(3) 



ksT 3 + exp{J/kBT) 

with a g-factor of 1.8, J = 1A6K and Xo = 8 x 
10~^emu/mol. We have re-analyzed the data by consid- 
ering the temperature dependence of the susceptibility 
for a spin-i alternating Heisenberg chain with: 



H — JlS2n • S2ri+1 + J2S271+I ' S2n+2)- 



(4) 



with Ji > and J2 > (antiferromagnetic cou- 
plings) obtained by the stochastic Quantum Monte Carlo 
methodEJ. We define the parameter a — J2/J1 which 
measures the ratio of the inter and intra-dimer coupling, 
a = corresponds to the dimer model. We introduce 
also an additional parameter 5 so that Ji — J(l -I- 6), 
J2 = J(l — S) and therefore a = 

In Fig. ^(a) we show the comparison of the experi- 
mental data to the dimer-model susceptibility given by 
Eq. H. We find, in fact, that the best fit corresponds 
to J = 1A6K, Xo = 8 X emu I mol and a value of 
g =r-J.77, somewhat smaller than the one proposed in 
Ref.Ea. Figs. ^ (b) and (c) show the comparison of the 
experimental data to the spin-^ alternating Heisenberg 
chain model for 5 = 0.8, a = 0.11 and 5 = 0.6, a = 0.25 
respectively. We observe that values of a up to 0.25 give 
still very good fits to the susceptibility data by choosing 
appropriate g and J parameters. This analysis makes 
plausible our ab initio results, namely the facL-jthat the 
inter-dimer interactions can be also significants. 

The important conclusion that we gain from this sus- 
ceptibility analysis is that it is a very hard task to deter- 
mine the behavior of a compound only from susceptibility 
fits. Extra information is needed so as to minimize the 
possible options. For instance, electron spin -resonance 
(ESR) experiments could complement the susceptibility 
measurements by delivering the value of the gyromag- 
netic factor. There are in the literature other examples 
{e.g. VOPO) where the only examination of the suscep- 
tibility led_tp a wrong interpretation of the behavior of 
the systemE3. 



We have presented ab initio band-structure calcula- 
tions for CSV2O5, a low-dimensional spin-gapped sys- 
tem and member of the alkali- vanadate family AV2O5 
(A=Na, Li, Cs) which has a similar monoclinic structure 
as the most studied VOPO. By means of a downfolding- 
tight-binding analysis, we observe that the behavior of 
this system can be best described by a spin-i alternating 
chain model with weak interchain interactions. This pro- 
posal has been complemented by an analysis of available 
susceptibility data for this compound. A re-examination 
of the temperature dependence of the experimental sus- 
ceptibility leads us to the conclusion that an alternating 
chain model can also provide a good fitting by choosing 
appropriate g and J values. We have discussed the struc- 
tural, electronic and magnetic similarities and differences 
between CSV2O5 and a'-NaV205, 7-LiV205 and VOPO. 

With this work we would like to draw the attention 
on this compound which shows interesting features, es- 
pecially since it has a lot of points of comparison to the 
most studied VOPO. We think that new experimental 
studies like ESR, INS, Raman scattering and the study 
of a possible field-induced magnetic ordering under a very 
strong magnetic field would be helpful to unveil the be- 
havior of this material. 
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FIG. 1. Crystal structure of CSV2O5 projected in the 
{yz) plane (see text for explanation). The large circles are 
the V-ions, black and grey for V*^ and V^^ respectively. 
The oxygens are represented by the smaller circles. Pairs 
of edge-shared distorted V^^Os square pyramids are bridged 
by V^'*'04 tetrahedra to form layers. The alkali-ion Cs shown 
by white circles , are located in between these sheets. 
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FIG. 2. DFT results for CSV2O5. The path is along 
r=(0,0,0), B=(-^,0,0), D=(-^,0,^), Z=(0,0,7r), F, Y=(0,7r,0), 
E=(-7r,7r,7r). Also shown in rectangles is the band character 
(in the local coordinate system) where V1=V^^ and V2=V^^. 
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FIG. 4. Hopping parameters corresponding to the 
tight-binding modehng of CSV2O5. Not shown is the on-site 
energy eo of the V^*''' site. While ti corresponds to a direct 
V^'^-O-V*"'' superexchange path, t2, t-i and are defined by 
pathes through the neighboring V^^04 tetrahedra. 
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FIG. 5. Comparison of the tight-binding bands (solid 
lines) with the DFT bands (dashed lines). The tight-binding 
parameters (see Fig. ^ are (in eV) eo = —0.0215, ti = 0.117, 
t2 = 0.015, t3 = 0.097, is = 0.050. 
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FIG. 3. Partial density of states of the V*^-3d orbitals 
in (a) CSV2O5 (b) i-LiVaOs (c) a'-NaVaOs obtained from 
DFT. Compare the distinct behavior at band edges around 
the Fermi level. 
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FIG. 6. Temperature dependence of the magnetic suscep- 
tibility of CSV2O5. rljhe dots correspond to the experimental 
data by Isobe et alxj where the Curie contribution has been 
subtracted. The solid lines show a fit to (a) spin-i Heisen- 
berg dimer with J/ks = 146 K, g — 1.77, (b) spin-i alternat- 
ing Heisenberg chain with S = 0.8, Ji/k_B = 146 K, g — 1.79 
and (c) spin-i alternating Heisenberg chain with 5 = 0.6, 
Ji/ks = 146 K, 5 = 1.81. 
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